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We report meson spectra obtained by using valence overlap fermion propagators generated on 
a background of 2+1 flavor domain wall fermion gauge configurations on 16 3 x 32, 24 3 x 64 
and 32 3 x 64 lattices. We use many-to-all correlators with Z3 grid source and low eigenmode 
substitution which is efficient in reducing errors for the hadron correlators. The preliminary 
results on meson spectrum, oq correlators, and charmonium hyperfine splitting for three sea quark 
masses are reported here. 
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1. Introduction 

During the last few years RBC and UKQCD collaborations have generated a large set of con- 
figurations with 2+1 flavor dynamical domain wall fermions (DWF) on several lattices with pion 
mass as low as 300 MeV and volume large enough for mesons {m n L > 4) [1]. Reducing sea quark 
masses to lower values requires substantial computational effort. A possible expedient approach 
toward unquenched QCD simulation with chiral fermions could be to use overlap fermions as the 
valence quark on those domain wall gauge configurations. In view of the fact that the overlap 
fermion has additional desirable features, such as incorporation of multi-mass algorithms and pre- 
cise evaluation of the matrix sign function, which one can take advantage of in order to improve 
chiral symmetry as well as the quality of the numerical results. 

The use of overlap valence on domain wall fermion gauge configurations is a mixed action ap- 
proach which is associated with both mixed action and partially quenched disadvantages. However, 
the mixed action partially quenched chiral perturbation theory (MAPQ^PT) has been developed 
for Ginsparg-Wilson fermion on Wilson sea [2] and staggered sea [3], and has been worked out for 
many hadronic quantities to next-to-leading order (NLO), such as pseudoscalar masses and decay 
constants [2, 3], isovector scalar oq correlator [4, 5, 6], heavy-light decay constants [7], and baryon 
masses [8, 6]. Using similar MAPQ^PT for two different chiral fermions, it is possible to extract 
reliable physical quantities such as masses and decay constants. Since both fermions used here 
are chiral fermions, A,„i x , the low energy constant representing G(a ) discretization dependence, 
is smaller than other mixed action formulations [10]. We shall also adopt the deflation method for 
overlap formulation, as will be mentioned later, to speed up inversion. The overlap formulation 
also incorporates multi-mass algorithm which enables calculation of multiple quark propagators 
covering the range from very light quarks to the charm on these sets of DWF lattices. This makes it 
possible to include the charm quark for calculations of charmonium and charmed-light mesons us- 
ing the same fermion formulation for the charm and light quarks [9]. With this formalism, while it 
is possible to get insight of light quark behaviors of many physical quantities, rich phenomenology 
involving charm mesons and baryons can also be studied in the same lattice formulation. 

2. Formalism : Overlap with deflation, HYP smearing, Z3 grid source and low mode 
substitution 

The formalism for our calculation was detailed in ref. [10]. Here we outline that briefly. It has 
been shown [11] that the projection of low eigenmodes from the Dirac operator can speed up inver- 
sion of fermion matrices and this procedure of deflation has been applied to both hermitian [12, 13] 
and non-hermitian [14] systems as well as to hermitian system with multiple right-hand sides [15]. 
In addition to faster inversion one can also substitute exact low eigenmodes in the noise estimation 
such as in quark loops [16] and all-to-all correlators [17, 18] to reduce errors in two and three-point 
correlators of mesons. 

One can get the eigenvectors of the massive overlap Dirac operator by calculating the same 
for the massless Dirac operator D ov . Since Dj v = ysD ov ys, and [D uv D\ v ,y^\ = one first uses 
Arnoldi algorithm to search for eigenmodes of D ov Dl v with real eigenvalues | A, | 2 which are doubly 
degenerate with opposite chirality. Eigenmodes of D ov can then be obtained by diagonalizing the 
two chiral modes in D ov . 
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Figure 1: (left (right)) The lowest 200 (400) eigenvalues for kernel in the inner (outer) loop of the overlap 
fermion for a 32 3 x 64 configuration with m; = 0.004. 
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Figure 2: (left) The vector correlator for the 32 3 x 64 lattice at m s <~ 200 MeV for three different choices, 
(right) The respective relative errors. 

Smearing the gauge links suitably could deplete the density of the lowest eigenvalues in Hw 
and by using HYP smearing we observe that the lowest eigenvalue with HYP smearing after defla- 
tion with 100 to 200 eigenmodes is about 3 times larger than those without HYP smearing. In Fig. 1 
we have plotted the spectra of the lowest 200 (400) eigenvalues for kernel in the inner (outer) loop 
of the overlap fermion for a 32 3 x 64 configuration with m/ = 0.004. In Table 1 we sum up our 
findings for speed up using HYP smearing and deflation in detail along with overhead of producing 
eigenmodes. 

Table 1: A comparison of speedup of inversion with HYP smearing (S) and deflation (D) of the outer loop. 
The inner and outer iteration numbers are for one propagator with 12 columns of color-spin. The speedup 
refers to that between the case of D + S vs the one with neither D nor S. 
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In addition to the above speed up procedure we also substitute the low-frequency part of the 
noise estimated correlator with the exact one obtained from the eigenmodes, which has been shown 
to reduce the variance [17, 18]. To do that we employ Z3 noise grid source with support on certain 
spatial grid points in a time slice to calculate the quark propagator which amounts to a many-to-all 
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Figure 3: The same as figure 2 for the nucleon correlator with m K ~ 300 MeV. 

approach as opposed to the all-to-all approach. This reduces the noise contamination in the high 
frequency part of the correlators. The details are given in ref. [10]. In Figs. 2 and 3 we have plotted 
the vector (at m n ~ 200 MeV) and nucleon (at m % ~ 300 MeV) correlators with point source, Z3 
grid source, and Z3 grid source with low mode substitution. On the right we also plot the relative 
eiTors which clearly show that Z3 grid source with low mode substitution is the better option. 

3. Numerical details 

We have calculated overlap quark propagators with overlap valence quarks by using 2 + 
1 flavor domain wall fermion gauge configurations generated by RBC and UKQCD collabora- 
tions [1]. Three sets of lattices are used with the four-dimensional sizes of 16 3 x 32, 24 3 x 64(a~ i = 
1 .73(3) GeV), and 32 3 x 64(a~ ! = 2.32(3) GeV) with several sea quark masses each. The Zolotarev 
rational polynomial approximation up to 14th degree is used to approximate the matrix sign func- 
tion and for the window [0.031, 2.5], the approximation to the sign function is better than 3.3 x 
10~ 10 [10]. Furthermore, low-mode deflation is used for both the inner and outer inversion of the 
overlap operator. Utilizing HYP smearing the largest absolute values of the overlap eigenvalues we 
deflate are 0.2, 0.125, and 0.22 on 16 3 x 32, 24 3 x 64, and 32 3 x 64 lattices with 100, 400, and 200 
eigenvectors, respectively. 

4. Results 

The results reported here are for the 32 3 x 64 lattice with lattice spacing aT x = 2.32(3) GeV. 
In Fig. 4(left) we plot {m n a) 2 as a function of valence quark masses ma in the range [0.0046 — 3.0] 
and for three sea quark masses. One can observe the nearly linear behavior up to a heavy quark 
mass. In Fig. 5(right) we also plot {m n a) 2 /ma vs ma for the same range of quark masses. One 
can observe the partially quenched chiral logarithm due to the mismatch of valence and sea quark 
masses. It is evident that the divergence becomes prominent with the increase of sea quark mass 
when the valence quark masses are lighter than the sea. In the future we will estimate the chiral log 
parameters for these sets of data. 

The effect of partial quenching is more evident in the case of the scalar meson. If the valence 
quark mass is lower than the sea mass then the lowest state in the scalar meson spectrum will be 
the would-be 77 n state which has negative contribution to the correlator. This has been investigated 
previously by various groups [21, 22, 5, 7, 6]. In Fig. 6(left) we plot scalar correlators for the sea 
mass 0.004 and for several valence quark masses. It is evident that as one increases valence mass, 
effect of the negative dip, which is a signature of the would-be r\n state, reduces. In Fig. 6(right) 
we plot scalar correlators for a given valence quark mass and for three sea quark masses. It is also 
evident that as one decreases the sea mass the contribution to the partially quenched ghost reduces. 
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Figure 4: (left) (m K a) 2 is plotted as a function of ma for the 32 3 x 64 lattice with quark masses ma = [0.0046 
- 0.3]. (right) (m n a) 2 /ma vs ma showing chiral logarithm for the same range of quark masses. 
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Figure 5: (left) Isovector scalar meson correlators are plotted for various valence quark masses and for 
m s = 0.004. (right)The same correlators at a fixed valence mass (lower than smallest m s ) and for various m s . 

If m vv i/m ss i is the mass of the pseudoscalar meson made up of the valence/sea quark and 
antiquark and m vs is the mass of the mixed valence and sea pseudoscalar meson, then the leading 
order (LO) pseudoscalar meson masses are given as 

m 2 w = B ov {m v + m v >), 
m 2 vs = B ov m v +B dw (m s + m TO ) + a 2 A mix , 

m 2 s , = B dw (m s + m s i + 2m res ) , (4.1) 

where A m , r is a low energy constant representing &(a 2 ) discretization dependence in the LO mixed- 
action chiral Lagrangian. A m « is also a measure of unitarity violation due to mixed action formalism 
at finite lattice spacing, which should vanish at the continuum limit. One can notice that A m i x 
term enters only in m vs . We have measured A miX term in ref. [10] and have found that it is ~ 7 
times smaller than the case of DWF valence on staggered sea [20] and ~ 1 8 times smaller than 
that of overlap on Wilson sea [19]. For a 300 MeV pion on the 24 3 x 64/32 3 x 64 lattice with 
a ~ 0.12/0.085 fm, the shift in mass due to A m j x is ~ 19/10 MeV. This small shift due to A m i x 
could admit a better extrapolation with the help of MAPQ^PT. 

In Fig. 7 we show meson masses obtained in a large range of valence quark masses at the 
lightest sea mass, m sea = 0.004. On the left side, we have plotted pseudoscalar, vector, axial and 
scalar masses as a function of {m n a) 2 . The threshold decay channels are also depicted with solid 
lines to indicate two-meson energy. For the vector channel, as one can observe, our data points 
are much below the threshold Tin P-wave energy indicating that Tin state lies higher than the p on 
on this lattice for the given range of quark masses. However, the axial channel decay threshold is 
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Figure 6: (left) Light quark meson masses are plotted over a range of (m K a) 2 . Threshold decay channels 
are also indicated by solid lines, (right) Meson masses in the charm quark mass ranges. 
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Figure 7: (left) The hyperfine splitting is plotted in the charm quark mass region as a function of -7= 
showing its linear dependence, (right) The same as a function of (m n a) 2 along with experimental data. 

overlapping our data which indicates that our data either represents the axial resonance {a\) or the 
two particle state p%. If this is indeed the pn state, then this will be the first observation of this 
state on the lattice. However, to conclude that one needs to show the existence of it along with 
the axial vector meson, and furthermore, the required volume dependence of a two particle state. 
On the right hand side we have plotted the meson masses in the charm quark mass range. In the 
future, we will do chiral extrapolation with MAPQ^PT and then continuum extrapolation to obtain 
physical values for these meson masses. 

In Fig. 8 we plot the hyperfine splitting (HFS) between heavy vector and pseudoscalar meson 
masses. Over the years, HFS obtained by lattice calculations are smaller than its experimental 
value. It was argued that lattice spacing errors along with dynamical fermion effects and non- 
inclusion of the disconnected insertion contributed to this discrepancy. The HFS is expected to scale 
like to leading order in m. Based on this expectation, in Fig.8(left) we plot the HFS obtained 
in the charm quark mass range for two sea quark masses. It is to be noted that the data points 
are linear in this mass range and one can expect a smaller m 2 dependence. We fit the HFS in this 
range to the form AEjjfs = 7/^(1 + wmc h includes the next term in large m expansion. The fit 
value obtained are a = 0.0349(5) and b = 0.0014(2268) which shows that the 1 / y/m dependence is 
pretty accurate. In Fig. 8(right) we plot the HFS as a function of m|a 2 along with the experimental 
number and one can observe that the experimental number is very close to one of our data point. 
Of course one needs to do continuum and chiral extrapolation to obtain its final value. This is very 
encouraging in the sense that other charm physics can also be obtained from this set of lattice. 
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5. Summary and outlooks 

Overlap valence on 2 + 1 flavor DWF configurations is a promising approach to do lattice 
QCD simulation with light, strange, and charm quark together in the same fermion lattice formula- 
tion. We have demonstrated that the eigenvalue deflation along with HYP smearing is very efficient 
procedure for inversion and the the Z3 grid source with low-mode substitution can reduce error in 
two-point correlator up to a factor 4. Results shown for the meson spectrum are very encouraging. 
With appropriate MAPQ^PT we will next do the chiral extrapolation and then continuum extrap- 
olation. Furthermore, we found that the 1 / yjm behavior in the HFS extends to the charm quark 
range and the result for the HFS is quite encouraging. This motivates us to study charm-light, 
charm-strange and charm-charm spectrum and other related phenomenologies in the future. 
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